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Various fidelity measures can be defined between two quantum processes especially when at least 
one of them is non-unitary. In this paper we consider two such measures of state averaged process 
fidelity, put forward an efficient procedure to compute them and investigate their mathematical 
properties as well as their relationship. 
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I. INTRODUCTION 

The destruction of classical determinism in the quan- 
tum limit necessitates a probabilistic description of the 
quantum reality. The understanding of quantum states 
and their evolution in a Hilbert space is an attempt to 
reconstruct the deterministic nature of quantum mechan- 
ics from a different perspective. Hilbert space in quan- 
tum mechanics plays the role of phase space in clas- 
sical mechanics and therefore the question of closeness 
between two quantum states becomes legitimate [H-[3]- 
However, for situations where quantum operations are of 
central importance, such as in quantum computing [l|, 
an equally relevant question is how to quantify the close- 
ness between two quantum operations. The requirement 
of defining a measure of fidelity between two quantum 
processes has so far been considered mostly for cases 
where the relevant quantum operations are unitary 0- 
.6], which motivates us to investigate a more practical 
question of introducing fidelity measures in the context 
of non- unitary quantum operations Q. 

The non-unitary quantum operations are relevant in 
various areas especially in superconducting quantum 
computing where the quantum gates are implemented 
within a subspace of the full Hilbert space [8[. In such 
cases the leakage outside the computational subspace 
constitutes an integral part of the gate error. The goal 
here is to introduce fidelity measures that can quantify 
this error efficiently. In this paper we consider two such 
measures. The first one is a natural generalization of the 
unitary case Q, where we project the full time-evolution 
operator (U /) first to the computational subspace (to ob- 
tain U s ) and then apply the fidelity definition for uni- 
tary operations between U s and our target quantum gate 
l/tgt- In the second measure we embed our unitary tar- 
get t/tgt inside a bigger unitary operator Vtgt having the 
same dimension as Uf and then maximize the fidelity 
(obtained via fidelity definition for unitary case) over the 
non-computational subspace of Vtgt ■ The motivation be- 
hind the second definition of fidelity measure stems from 
the fact that in a real experiment we produce the en- 
tire time evolution operator Uf, but measuring the error 



only in the computational subspace tacitly assumes that 
the error in the non-computational subspace is minimum 
if not zero. This can be achieved only if the fidelity is 
maximized over various choices of constructing a target 
in the non-computational subspace of Vt g t ■ In this paper 
wc explore the mathematical properties of these two fi- 
delity measures and derive a relationship between them. 
We show that the expensive computation involved in the 
maximization of the second definition can be overcome 
by obtaining the maximum value analytically applying 
the polar decomposition theorem. Our numerical result 
reveals quantitative features of these fidelity measures for 
a realistic situation. 

The rest of the paper is organized as follows. In Sec. pl)) 
we first describe the fidelity measure for unitary opera- 
tions and then introduce two measures of fidelity of non- 
unitary quantum operations. In Sec. we analyse var- 
ious mathematical properties of the fidelity measures. In 
Sec. (|iV[) we perform an analytical comparison between 
the two definitions. In Sec. fVT) we discuss how our re- 
sults can be generalized for cases where the target quan- 
tum operation is also non-unitary. We demonstrate our 
analytics with a numerical example in Sec. (|VI|) and dis- 
cuss the conclusion and scopes of our work in Sec. dVIII) . 



II. FIDELITY MEASURES 

In this section, we discuss fidelity measures for unitary 
and non-unitary time-evolution operators. 



A. Unitary Operation 

Fidelity between two quantum states denotes their dis- 
tance in the Hilbert space and can be defined as, 



F(\ip),\ X )) = \mx)\' 



(1) 



where and \x) are two (pure) states. We can always 
define a state dependent fidelity between two quantum 
processes (U and Ut) as, 
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F X (U U U)= (X\U}U\ X ) , 



(2) 
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where \x) is an arbitrary quantum state on which the 
quantum operations are applied. For cases when U and 
Ut are both unitary, the state averaged process fidelity 
between them can be defined as, 



F avs (U u U) = 



(x\u}u\ x ) 



dV, 



(3) 



where \x) is a (n dimensional) normalized state vector 
defined on the unit sphere S 12 ™ -1 in C and dV being a 
normalized measure. It can be shown 0, @] that such a 
state averaged fidelity formula can also be expressed as 
a trace formula given by, 



F avg (U t ,U) = 



Tr(U}uWU t )+ Tr{U}U) 



n(n + 1) 



(4) 



n being the dimensionality of the Hilbert space. While 
the result in Eq.Q being a mathematical identity irre- 
spective of the unitarity of U and Ut, the definition in 
Eq.Q is necessarily applicable for situations when both 
Ut and U are unitary. Such cases arise when we compute 
the fidelity between two quantum operations defined on 
a full Hilbert space or a subspace which is absolutely de- 
coupled from the rest of the Hilbert space. In quantum 
computing such cases are highly unlikely and therefore 
the interesting question is how to define a fidelity mea- 
sure between two quantum operations where the actual 
time-evolution operator U may be non-unitary while our 
target quantum operation Ut g t is unitary. 



B. Non-unitary Operation 

There exists at least two different ways to define a fi- 
delity measure for non-unitary cases and in this paper 
we investigate their mathematical properties and the re- 
lationship between them. In quantum computing usually 
the non-unitary cases arise when we have auxiliary states 
3 or coupled subsystems We are interested to com- 
pute the fidelity of a quantum operation (U s ) which lies 
in a subspace of the full Hilbert space and therefore non- 
unitary in general. If we denote the full time-evolution 
operator (unitary) defined on the entire Hilbert space by 
Uf, then in some matrix representation, 



Uf 



B | C 

V I 



A \ 
I 



(5) 



where U s is m dimensional, Uf is n dimensional and A, 
B, C are other blocks of the full time-evolution operator 
(Uf). The only constraint we require here is Uf to be 
unitary. 

One way to define fidelity is to construct U s from the 
full time-evolution operator Uf projecting it to the rel- 
evant computational subspace and then compute the fi- 
delity between U s and the target quantum operation Ut g t 



as it is defined in Eq.Q, 
F^Ut^Us)^ 



(x\ uLu, \ x ) 



Tr(U s Uj) + Tr(U*U B ) 



dV 

2 



m(m + 1) 



(6) 



where both U s and Z7 tgt are m dimensional as shown in 
Eq.©. 

Another way to define fidelity is to embed the target 
quantum operation Utgt (m dimensional) inside an n di- 
mensional unitary matrix Vt g t as, 



Vi 



t g t 



V 



K 



3 \ 
) 



L 



(7) 



and then choose J, K, L such that the fidelity (as de- 
fined in Eq.©) between Uf and Vtgt is maximum. This 
maximisation is required because in a real experiment 
the error of the gate operation performed in the non- 
computational subspace is implicitly assumed to be min- 
imum and therefore it gives an alternative fidelity mea- 
sure between Utgt and U s and is given by, 



Tr{U f U) 



F 2 (U tgU U s ) = F^U tgt7 U s )\ 
F^(U tgU U s ) = 

+ Tr(VhUff 



max{ J,K,L} •> 
2 



(x\ v t f gt u f \ x ) 

Tr(V^U f ) 



dV 

2 



n(n + 1) 



n(n + 1) 



"(8) 



where max{J,K,L} means we choose a set {J,K,L} for 
which F 2 is maximum. The last equality holds true as 
Uf is always unitary. 



III. ANALYSIS 

First, we explore the mathematical properties of these 
definitions. It's interesting to note that if we require 
both Utgt and Vtgt to be unitary then it's possible only 
if J = K = and L is unitary. This enables us to write 
Vt g t as a direct sum of Utgt and L as, 



Vtgt = Utgt © L 



(9) 



We can write the 



Tr{Vt\tU f ) 

Tr(V t lU f )\ = \Tr{U^U.) + Tr(L^C) . (10) 
Let's denote, 

Tr(uLu s ) ~ re w and Tr(tfC) ~ se 4 * (11) 
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Now, we can write F!± in terms of r, s and F\ as, 



+ m(m + l)Fi - Tr{U s Ul) + s 2 + 2rs COS (V — 
n(n + 1) 



(12) 

where s = \Tr(lJC)\. We show that |Tr(L t C)| must 
have a finite upper limit using the following theorem. 

Theorem 1: If X (not necessarily unitary) be a 
pxq dimensional block inside a unitary matrix W , then 
Tr{XX^) ^ mm{p,q}. 

Proof: Let's write the unitary matrix W as, 



X ! Y 



W EE 



V 



\ 



/ 



Since W is unitary, Tr(XX^) + Tr(YY^) = p and also 
Tr(X^X) + Tr(Z^Z) = q. For any matrix T, Tr(TT t ) 
is a positive real number. Therefore, Tr(XX') ^ p and 
also Tr{XX^) < q, that altogether gives Tr(XX^) ^ 
min{p, q}. Physically the maximal case occurs when the 
entire population leaks from outside to the blocked sub- 
space. 

For our purpose since L is unitary, Tr(L^ L) = n — m 
and using Theorem 1, Tr{C*C) ^ n — m and therefore 
(using Cauchy-Schwarz inequality) we can show, 



s ee \Tr(UC)\ jTr{L^L)Tr{C^C) < »- 



(13) 



It is also possible to determine the maximum value of s 
(smax) by using the polar decomposition theorem. Ac- 
cording to polar decomposition theorem any square ma- 
trix can be written as a product of a unitary and a 
hermitian matrix, where the hermitian factor is always 

I 



unique. So we can write [10( C as, C = CtjCh where 
Ch = VC^C, the square root here denoting the principal 
square root of a matrix where all the eigenvalues have a 
non-negative real part. Since the Cjj is the closest uni- 
tary matrix to C [lOj . we can choose our L to be equal 
to Cu and write, 



Tr(ClC v Cn) = Tr(VcW) , (14) 



where the square root denotes the principal square root 
of a matrix. 

Now as we can see from Eq. (I12p , in order to maximize 
F' 2 we need to satisfy two conditions: i) we need to choose 
an L for which s is maximum (given by Eq. (|14p ) and ii) 
we need to choose an L for which = <f>. Suppose L max 
satisfies these conditions and s max being the maximum 
value of s. We can define a fidelity between C and L as, 



Font (L,C) = F avg (L,C) 



Tr(C^C) + \Tr(tfC)\ 
(n — m)(n — m + 1) 



(15) 

We can see from Eq.(fT"5j) that s is maximum iff F out (L, C) 
is maximum, but we can't find a unitary L for which 
-Fout is unity if C is non-unitary. In such situations i max 
actually denotes a unitary matrix closest to non-unitary 
matrix C (which is Ctj as defined earlier). If F™ x is the 
maximum possible fidelity between L and C, then from 
Eg . (TTTl and [TBI . we can write, 



Tr{C^C)+ Tr{VcU3) 



(n — m)(n — m + 1) 



(16) 



Now we write our second definition of fidelity in terms of 
first one as, 



Tr{VcW) +m(m + l)F 1 -Tr(U s W s ) 2 Tr(VCW) \/ m(m + l)Fi - Tr(U s u!) 



n(n + 1) 



n(n + 1) 



(17) 



which means F\ and F2 carries a non-linear relationship 
between them. 



IV. IMPLICATION 



Eq. (fT7|) gives a neat formula between two definitions of 
fidelity. In this section we compare these two definitions. 
First, we show analytically that F2 F\ if U s is uni- 
tary. This actually means we need to prove the following 



inequality, 



Tr{Ul Rt U s ) 



Tr(uLu.) 



(18) 

since if U s is unitary then so is C. Notice that for this case 
(using Cauchy-Schwarz and Theorem 1), 



Tr{uyj. 



by 5 and therefore, 



m. Let's denote yn — Tr(U^ gt U s ) 

^ 5 < m. So, in order to show Eg. (ITS)), we need to 
prove the following inequality for any 5 ^ m, 



m + (m — <5) 2 n + (n — S) 



(19) 
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Let's start with the fact that n > m which implies 
11 ^ m + 1 and after a little manipulation this gives 

Iran 



to ^ , 

^ m+n+l 

that O^Km 



Combining this result with the condition 
we write, 6 ^ _ 2 ,'"" , . An algebraic ma- 



m+n+l ' 



2m-5 



nipulation over this step actually gives us m+m -i ^ 



2n-(5 



which ultimately gives rise to Eq. (|19[) and thus we can 
prove Fi ^ F\ when U s is unitary. 

Now in order for everything to make sense we need to 
show that at higher values of F\ there should not be much 
difference between Fi and F2 and we show here that at 
limj^^.i-Fa = 1. At Mmp^i, U s tends to be a unitary 
matrix and so is C and therefore plugging everything in 
Eq. flTT)) . we obtain, 



lim F2 



n + (n — m) 2 + m(m + 1) — to + 2(n — m)m 



n[n + 1) 



= 1. 



So, both the fidelity measures are identical at unity. 



V. NON-UNITARY TARGET 

It is possible to generalize our analysis for non-unitary 
target operations. Such a situation is unlikely in the 
context of quantum computing (without decoherence) , 
but may arise in the context of quantum simulation. By 
non-unitary target operation we actually mean our ZTtgt 
to be non-unitary while we still require Vtgt and Uf to 
be always unitary and as we can note for such cases J 7^ 
K 7^ in general. This essentially means that we need 
to replace L^C by J^A+ K^B + L^C in our previous 
analysis which gives, 



+ m{m + l)Fi - Tr{U s Ut) + s 2 + 2rscos(0 - 



n(n + 1) 



where, 
Tr(U*U a ) 



(20) 



jo 



and TriJ^A + K^B + L^C) 



(21) 

for this case. In order for F' 2 to be maximum we need 
to choose J, K and L such that s is maximum. So, it's 
important to show that such a maxima exists. Applying 
triangle inequality we write, 



s < \Tr(,PA)\ + \Tr(K*B)\ + \Tr(L f C)\ 



(22) 



And finally using Cauchy-Schwarz inequality and our 
Theorem 1, we can show that s ^ m + n; m, n being 
dimensions of C/ tg t and Vtgt respectively. This actually 
means that s must have a finite upper bound and sup- 
pose s max be its possible maximum value. Then we can 
write, 



Fi 



■ m{m + l)Fi — Tr(U s Ul) + s, 2 nax + 2rs 
n(n + 1) 



(23) 

which is the relation between two fidelity measures for 
non-unitary target operations. 



VI. NUMERICAL ANALYSIS 

In this section we investigate the two measures of fi- 
delity more quantitatively in the light of a real physical 
system. We consider a controlled-Z gate implementation 
in a superconducting phase qubit capacitively coupled to 
a resonator. The Hamiltonian for this system is given by, 



H = Hq + H u 



(24) 



where, 









UJ r 


) 


2w r 








/ -i 





I i - 




V y/2i 






( 25 ) 

where, cj q , uj v and A are qubit frequency, resonator fre- 
quency, and anharmonicity of the qubit, respectively and 
g is the (time-independent) interaction strength between 
the qubit and resonator. The Hamiltonian used here has 
nine levels in it while four of those constitute a computa- 
tional subspace. So we can define and compute both the 
fidelities for this situation. The purpose here is to com- 
pare two fidelity measures quantitatively and in order to 
simplify our problem we set the qubit frequency w q (6.2 
GHz.) at uj r + A, resonator frequency uj r at 6 GHz. and 
coupling (g) at 30 MHz. and consider our Hamiltonian 
to be time-independent. 




9.76 9.78 9.8 

t (in nanoseconds) 



9.82 



FIG. 1: (Color online) A plot of two fidelity measures 
vs. time (in nanoseconds) where the target operation is 
a controlled-Z gate in a resonator coupled phase qubit. 

We consider the time evolution of our Hamiltonian for 
a short span near the regime where we expect to obtain a 



5 



controlled-Z gate (Strauch's approach and plot both 
the fidelities (Fig. (JT])) with respect to time. Our numer- 
ical result shows that F2 > F% for that regime. We have 
proved analytically that F% > F± if U s is unitary, while 
proving (or disproving) the result Fi F\ universally 
(for non-unitary situation) has not been accomplished in 
our work and therefore is left as an open problem. 



VII. CONCLUSION 

In summary, we have considered two possible measures 
of fidelity for non-unitary quantum operations and ex- 
plored their mathematical properties. A relation between 
them has been established. It is found that the maxi- 
mization involved in the second definition of fidelity can 
be done analytically and the optimal value is expressed 
in closed form. Our numerical result shows the quantita- 
tive feature of the two fidelities for a realistic situation. 
The plot in Fig.([T]) shows that while the two measures of 
fidelity differ considerably for lower values, they converge 
for higher values. The difference between two measures 
of fidelities at lower values is not surprising and shows 
up even for unitary case. For example, if we look into 
the relationship between the fidelity measures (for uni- 
tary situation) considered by Hill (let's call it Fh- Sec 
Eq.(13) in Ref.Q.) and the one considered by Pedersen 



et. al. (Let's call it Fp. See Ref. j3>]), we can write, 



F F 



nFl 



1 



(26) 



n being the dimension of the Hilbert space. Eq. (j2l)|) tells 
us that at F H 1, F P >-> 1, but at F H >-> 0, F P >-► j^. 
The difference between fidelity measures for lower fideli- 
ties resemble the fact that in statistics different measures 
of central tendency (mean, median and mode) differ a 
lot in general, while give identical values when all num- 
bers the data set are equal. Our numerical results also 
show that F2 ^ F\ . The proof (or contradiction) of this 
result for non-unitary cases is still lacking and can be 
considered as a future direction of research on this sub- 
ject. Our result only shows the relationship between two 
possible measures of state averaged process fidelity when 
a non-unitary time-evolution operator is produced in an 
experiment. Similar questions can be investigated as a 
follow-up of this work for state-minimized (worst case) 
process fidelity which is more relevant to quantum error 
correction. 
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